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Abstract. 

We present a covariant wave-packet approach to neutrino flavor transitions in 
vacuum. The approach is based on the technique of macroscopic Feynman diagrams 
describing the lepton number violating processes of production and absorption of 
virtual massive neutrinos at the macroscopically separated spacetime regions ( "source" 
and "detector" ) . Accordingly, the flavor transitions are a result of interference of the 
diagrams with neutrinos of different masses in the intermediate states. The statistically 
averaged probability of the process is representable as a multidimensional integral of the 
product of the factors which describe the differential flux density of massless neutrinos 
from the source, differential cross section of the neutrino interaction with the detector 
and a dimensionless factor responsible for the flavor transition. The conditions are 
analyzed under which the last factor can be treated as the flavor transition probability 
in the usual quantum mechanical sense. 
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1. Introduction 

The neutrino-oscillation hypothesis suggested by Bruno Pontecorvo pj offers the most 
probable and promising explanation of the anomalies discovered in many experiments 
with solar, atmospheric, reactor and accelerator neutrinos and antineutrinos. However 
the conventional quantum-mechanical (QM) theory of neutrino oscillations faces a 
number of methodological difficulties and paradoxes which have been widely discussed by 
many authors (see, e.g., [2HH] and numerous references therein). It has been recognized 
that, remaining within the framework of quantum mechanics, some problems of the 
conventional theory can be resolved by using wave packets to describe the neutrino 
states and/or the states of the particles involved in the processes of neutrino production 
and detection. In this way it has been found, e.g., that neutrino oscillations vanish if 
the neutrinos propagate over the distances much larger than the coherence length (for 
in-depth discussion, see |6]). 

A more radical approach to the problem is based on quantum field theory 
(QFT) P4T3]: see also [6H8] for further references. In this approach, the sequential 
processes of neutrino production and detection are treated as a single process which 
can be described by the S'-matrix formalism of QFT. The peculiarity (or better to say 
novelty) of this approach is in the use of Feynman diagrams with vertices separated by 
macroscopically large spacetime intervals. Besides, as in the QM case, for the description 
of the asymptotically free "in" and "out" states one has to use wave packets rather than 
the ordinary (in QFT) Fock states with definite 4-momenta. In this framework, the 
neutrinos appear as the propagators of the neutrino mass eigenfields, z/j, connecting the 
macroscopically distant parts of the Feynman diagram. In the new approach, there is no 
need in even mentioning the flavor neutrino states or fields, and the neutrino oscillation 
phenomenon is a reflection of the usual interference of the Feynman diagrams with 
different i/j's in the intermediate states. Such treatment not only enables us to reproduce 
the standard QM formula for the flavor transition probability and obtain the pertinent 
QFT corrections (thus offering a potential possibility for an experimental verification of 
the theory) but also impels us to re-examine the concept of neutrino flavor mixing. 

The authors of the mentioned pioneer papers do use (explicitly or implicitly) the 
noncovariant, usually Gaussian wave packets, which are only applicable, as it will 
be shown in the next section, to nonrelativistic particles. However, in the reactions 
of neutrino production and detection taking part in the current neutrino oscillation 
experiments, relativistic particles rule involved||] Moreover, by an appropriate 

Lorentz boost, any particle can be made to be ultrarelativistic. Therefore it is necessary 
to work out a relativistically covariant (inertial frame-independent) theory of wave 
packets applicable for the description of particles with arbitrary momenta. 

In this paper we propose a simple covariant wave-packet theory (section [2]), which is 
then applied to the calculation of the interference of macroscopic diagrams with massive 

I An important exception from the rule might be experiments with Mossbauer antineutrinos which are 
under active discussion in the current hterature p^UIS] , 
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neutrino exchange (section S]). We perform a statistical averaging of the probabihty (the 
squared absolute value of the overall amplitude) which leads to a rather general formula 
for the neutrino event rate in the detector (section [5l); the conditions of applicability of 
the conventional QM formula for the flavor transition probability are also discussed. In 
section |6] we summarise the main results of the formalism. 



2. Relativistic wave packets 

2.1. General consideration 

The ^'-matrix formalism of QFT usually deals with one-particle Fock states 



|fc,.) = V^aLlO) (1) 

as asymptotically free states of the fields. Here |0) is the Fock vacuum state (afcs|0) = 0), 
a^^ and a^s are the creation and annihilation operators of a particle with 3-momentum 
k and spin projection s, Ef, = + and m is the mass of the particle. The 
conventional (anti) commutation relations hold: 

where the curly brackets denote a commutator (for boson fields) or anticommutator 
(for fermion fields). The Lorentz- invariant normalization of the states ([1]) is therefore 
singular since 

(qr, r\k, s) = {2nf2Ek6srS {k - q) . 

Let us define the wave-packet state localized in a spacetime point x = {xq, x) as a linear 
superposition of the one-particle states ([1]) 



in which (j){k,p) is a Lorentz-invariant function. The last requirement implies that the 
proper Lorentz transformation A of the 4- vectors p = {po,p) and x = {xq, x), 

p I — > p' = Ap, X I — > x' = Ax, 

induces the following unitary transformation Ua of the wave-packet state: 

\p,s,x) I — > Ua\p,s,x) = \p',s,x'), 

assuming that the axis of spin quantization is oriented along the boost or rotation axis. 
We see that this transformation rule is in fact the same as that for the Fock states (see, 
e.g., [I9]) and this is just the reason of the form of definition given by ([2]). 

We will suppose that the function 0(fc, p) has a unique sharp maximum at k = p, 
and its behaviour around the maximum is governed by a small parameter a of the 
dimension of massl§| We further require that the wave-packet state ([2]) passes into the 

§ In general, the function (j){k,p) may involve a finite or infinite set of parameters. Here, solely 
for simplicity, we assume that [dln0(fc,p)/(i(fc — cx cr~^ > 0. Then a is the only essential 

combination of the mentioned parameters. 
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Fock state ([T]) as a — > 0. This is possible if 

lim (p{k,p) = {2tx f2Ep5{k - p) . (3) 

tr— >0 

This correspondence principle enables us to impose on the function 0(A;, p) the following 
Lorenz-invariant condition: 

dk(j){k,p) f dk(j){k,0) 



{2T:f2Eu J (27r)32Efe 
Indeed, the dimensionless Lorenz-invariant integral (jlj) does not depend on p. Hence it 
can only depend on the dimensionless ratio a/m. Owing to the integral (jl]) tends 
to 1 as cr — 0; thus it is natural to put it equal to one also at finite (small) a. Finally 
it can be said that the form factor function (j){k,p) represents, up to a multiplier, a 
"smeared" 5 function in the momentum space. 

Now, we find the wavefunction describing the state ([2]) in the configuration space. 
Let us do this for a spin- 1/2 fermion with the field operator 



The state (0|\l/(x) can be treated as a linear superposition of one-particle states with 
definite momenta at the point x which cannot be characterized by a specific momentum. 
The coordinate representation of the Fock state ([1]) is given by its projection onto the 
state (0|^(a;), 

(0|*(a;)|p,s)=M,(p)e-^P^ (5) 

which is a plane wave uniformly distributed over the whole spacetime. In contrast, the 
wave-packet state (|2]) is characterized by a momentum distribution governed by the 
function (j){k,p) concentrated near the most probable momentum p. Besides that, the 
spacetime distribution of the wave packet is not uniform as can be seen from the analog 
of equation ([5]) for the state ([2]): 

(0|^(x)|p, s,y)= e-'Py [u,ip) - VpM,(p) ■ («V. + p) + . . .]^ip,y - x), 

^ e~'''Pyus{p)ip{p, y-x). (6) 

Here we have introduced the Lorentz-invariant function 

/dk f dk 

which satisfies the Klein-Gordon equation and describes the wave packet in the 
coordinate representation. The last equality in ([7]) is written in the rest frame of the 
packet {p^ = 0); the 4- vector = (x^, a;^) is related to x = (x*^, x) through the Lorentz 
boost along the direction —p. Since ip{Q, x^,) is an even function of a;^, it can depend 
only on the variables x° and Icc^l, related to the invariant quantities (px) = mx^ and 
x^ = xl- The approximation ([6]) is valid under the condition 

\iWx\Y).ip{p,y - x) +p| < 2Ep, 

which is fully consistent with other approximations in the subsequent analysis. Note 
that the approximate equality (jH]) becomes exact for (pseudo) scalar fields. In fact, with 
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a certain loss of simplicity and universality of the formalism, it would be possible to 
completely get rid of this approximation if dealing with a spinor wavefunction instead 
of the scalar function ([7]). 

From ([7]) and ([3]) it follows that ip{p,x) — )■ e*^^ as cr —)■ 0, i.e., as expected, a wave 
packet utterly localized in the momentum space is delocalized (spread-out) over the 
configuration space. The inner product of the states ([2]) is nonsingular and is given by 

(q, r, y\p, s, x) = 5,,e^(")p(p, q;x-y), (8) 

where 

r rlU 

V{p.q;x) = J ^^^^3^^^ 0(fc,p)0-(fc,q)e-^-. (9) 

From dH]) and (Q it follows that the normalization of the state ([2]) is also finite: 

{p,s,x\p,s,x)=Vip,p;0) = 2EpV{p). (10) 

The quantities Ep and V(p) in (ITU!) are, respectively, the mean energy and effective 
spatial volume of the packet, defined by 

— _ J dxi:{p,x)idoij*{p,x) _ 1 C dk\(j){k,p)\'^ 

Jdxmp,x)\^ V{p)J A{2nrEk ' ^ ^ 

where Tp = Ep/m. So, both Ep and V(p), as well as the mean momentum p defined 
by a relation similar to f lTT]) . are integrals of motion. It is easy to prove that the mean 
position of the packet follows the classical trajectory: 

X= j dxi)*{p,x)xi){p,x) = VpXo. (13) 

Here Vp = p/Ep is the mean group velocity of the packet, which coincides with the 
most probable velocity VpEp = p/ Ep. 

Certain properties of the function (jH]) become especially transparent in the center- 
of-inertia frame of the two packets (defined by the condition + 9=,, = and denoted 
from here on by an asterisk subscript). In this frame 

/dk 
— — — 0(fc, pjr (fc, -pje^^(--^-). (14) 
[z-KyzEk 

Due to the assumed sharp maximum of (j){k,p) at k = p, one may expect that the 
function (1141) has a sharp maximum at = (that is at g = p) and quickly vanishes 
at large values of |p^| since the maxima of the multiplicands \(f){k,p^)\ and \4>{k, —p^,)\ 
in the integrand are widely separated in this case and thus their product is small for 
any k, as is illustrated in figure [H One can see in addition that the integral (fT4ll 
vanishes at any if the points a;* and are sufficiently separated in space (namely, 
if \x^: — ^ 1/^) since the phase factor Q-^H^*-y*) in the integrand of ([HI) rapidly 
oscillates in this case. 
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Figure 1. Schematic illustration of the vanishing of the integrand in equation at 
large 



2.2. Relativistic Gaussian packets 

In further consideration we will use a simple model of the state 
packet (RGP), for which the function (f){k,p) is of the form 

27r2 / EkEp-kp 



relativistic Gaussian 



(j){k,p) 



exp 



dcf 



'k,p), 



(15) 



where Ki{t) is the modified Bessel function of the third kind of order 1. One can verify 
that the function (ITSl) has the correct plane- wave limit ([3]) and satisfies the condition 
In what follows we assume ^ m^. Then the function (fTSl) can be rewritten as an 
asymptotic expansion by a small parameter a'^/m'^: 



[k,p) 



m 



a" 



a 



exp 



1 + 



3a' 
4m2 



(16) 



It is easy to convince that the nonrelativistic limit of the function f lT6|) coincides, 
up to a normalization factor, with the usual (noncovariant) Gaussian distribution 
<^g(^ ~ p) oc exp [— (fc — p)^/4cr^] widely used in the literature. However at relativistic 
momenta the functions 0g and (fc significantly differ from each other. For example, in 
the vicinity of the maximum k = p, 



•^G 



{k,p) 



27r3/2 



m 



a 



exp 



[k 



P) 



[k P). 



We see that in this case the relativistic effect consists in a "renormalization" of the 
wave-packet width {a — > aVp). This renormalization is very essential for the neutrino 
production and detection processes involving relativistic particles. 

The coordinate representation of the RGP wavefunction is found to be 



^g{p,x) 



Ciri(mV2o-2) 



(17) 
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where ( is the complex dimensionless scalar variable: 

C = W 1 [a'^x^ + i{j)x)\. 

V 

Here and below the square root means the principal square root. Figure [2] shows the 
shape of the function |'(/'g(0, a;^)| calculated, as an example, for ajm = 0.1 and plotted 
as a function of the dimensionless variables a'^x^/m and a'^xl/m, assuming that x^, is 
directed along the third axis; of course, \iPg{0,X-),) \ is an even function of both variables. 
As one can observe, |'?/'g(0, Xi,) \ rapidly vanishes with increase of x^ and its spatial width 




Figure 2. A 3D plot of 1-00(0, a:*)| as a function of a^x^/m and a^x^/m (assuming 
tliat Xi, = (0,0, a;^)). Tlie calculations are done for a/m = 0.1. 

grows with x^, i.e. RGP spreads with time as any wave packet with nonzero mass. 
Nonetheless the effective volume (fT2l) does not depend on time since the spreading of 
|'?/'g'(0, a;^)| precisely compensates its fall-off and the normalization of the RGP state 
remains constant. 

In the RGP model, the function ([9]) which defines the scalar product of the states 
dH]) reads 

q- x) = iE^-^ip fl^T'Jfl = q; x); (18) 

the complex dimensionless scalar variable z in (fT8|) is defined by 

z = - — \/ (p + qY — 4cr2 [cr^a;^ + iin + q)x\. 
2m 

Now we have to determine the physical conditions under which the spreading of 
RGP can be neglected, since just in this regime the wave packets can naturally be 
associated with the (quasi)stable particles and then used (instead of the plane waves) as 
the asymptotically free states of incoming and outgoing fields in the S'-matrix formalism. 
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2.3. Nondiffluent regime 

An accurate analysis of the asymptotic expansion of In [V'gIO, x^)] in powers of the small 
parameter a^/^ni^O, taking into account the inequalities \(\ > 1 and |arg(^| < 7r/2, 
provides the following (necessary and sufficient) conditions of the nondiffiuent behavior: 

a\xiy « mVc^^ a^\x,\^ « m^a^ (19) 
They can be rewritten in the equivalent but explicitly Lorentz-invariant form: 

(px)^ < m7(x^ [px)^ - m^x^ < m^ja^. (20) 
Under these conditions, the function (fT7|l reduces to the simple and transparent form 



(21) 



V'g(P! a;) = exp {%mx\ — a^a^^) = exp |i(pa;) ^ \^V^)^ ~^ m^x^] | . 

Let us mention some properties of this approximation which will be referred to as the 
contracted relativistic Gaussian packet (CRGP). It is easy to verify that the mean 
coordinate of the packet follows the classical trajectory (fT3l) and the absolute value of the 
function ipcip, x) is invariant under the transformations xq i — > xq + r, x i — y x + VpT. 
It is also obvious that \iPg{p,x)\ = 1 along the classical trajectory x = VpXo and 
|'?/'g(p, x)| < 1 with any deviation from it. In the nonrelativistic limit, the wavefunction 
(EH) takes the form 



i/^g{p, ~ exp [im {xq — Vpx) — a-"^ \x — VpXof'\ . 

In the CRGP model one obtains that V(0) = [n / {2a^)]'^/^ = and Ep = Ep. 

The CRGP approximation for the function f[T5]) can be derived by analyzing the 
asymptotic expansion of its logarithm in powers of the parameter o'^ j{vp?z) and taking 
into account the inequalities 1^1 > 1 and | arg z\ < 7r/2. This yields 



1^g{p*, -P*; X,) = exp 



(r* - 1) cr^x: 



,2^2 



(22) 



a2 2r, 

where F^, = E^^/m [E^, = Ep*). The applicability conditions of this approximation, 
o-^(x°)^ ^ El/(y'^ and cr^|a;*p ^ Eljo'^, are found to be fully compatible with 
the conditions (fT9|) . As was expected (see section I2.2p . the function (122|) rapidly 
vanishes if either or \x^\ (or both) are sufficiently large. The function Dq has 
a number of unobvious a priori properties. For example, exponentially vanishes at 
subrelativistic energies (F^, — 1 ~ 1). At nonrelativistic energies, the function Dq in the 
lab. frame is given by 



Vcip, q; x) ^ 2mV^ exp 



im {xq — vx) 



m'^\Vp — Vq\ a'^ \x — vxo\ 



8(t2 



where Vp = p/m, Vq = q/m (assuming that \Vp^q\ <^ 1) and v = j {Vp + Vg). 

As an application of equation ( l22l) . we consider the norm of a state with two identical 
noninteracting packets. The following exact model-independent relation holds: 

{Pl,Si,Xi;P2,S2,X2\Pi,Si,Xi;P2,S2,X2) . , ^ \V{p^,P2;Xi - X2)\'^ 



[2EpVipW [2EpY{pW 
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■'— " refer to bosons and fermions, respectively. From 
= ^2 = 0' ^he right-hand side of fl23|) is equal to 
. Hence the effects of Bose-Einstein attraction and Pauli 
S2) are essential only if the distance between the centers of 



where the signs "+" and 
f l22|) it follows that at 
1 ± Ss^s2 exp i-a'^\xi - X2\^] 
repulsion (appearing for si = 
the packets is comparable or less than their effective dimensions {\xi — X2\'^ ^ 



that 



is, exactly when it is already necessary to take into account the dynamic interactions 
between the packets. At sufficiently large distances between the packets, the quantum 
statistics is of no importance any more. Similar conclusions can be proved to be valid 
for the states with arbitrary number of free identical packets. 

3. Overlap integrals 

We will deal with the Feynman diagrams, the generic structure of which is shown in 
figure [3l The external legs of such diagrams correspond to asymptotically free incoming 



X 



y 



X, 



Figure 3. A generic macroscopic Feynman diagram with exchange of massive neutrino. 



("in") and outgoing ("out") wave packets in the coordinate representation, that is, to 
the wavefunctions ipa iPai^a) and ipl {pfj,Xb) specified by the most probable momenta 
p„ ft, spacetime coordinates Xa,b, masses m^^ft and parameters aa,b- Here and below we use 
the following notation: Is (Fg) is the set of in (out) packets in the block Xg ("source"). 
Id {Fd) is the set of in (out) packets in the block X^, ("detector"). The internal line 
connecting the blocks Xg and X^, denotes the causal Green's function of the neutrino 
mass eigenfield {i = 1,2,3). The blocks Xg and Xd are assumed to be macroscopically 
separated in spacetime; this explains the term "macroscopic diagram". 

In the calculations with the macrodiagrams of such kind, we will encounter the 
four-dimensional overlap integrals Yg{q) and Vrf(g) defined as follows: 



V.,,(g) 



dxe 



±iqx 



X] 



n 



e^^^^V,* (Pft, Xb 



In the CRGP approximation, these integrals can be written in the form 



x] 



VM(g) 



dxexp 



I {±qx - qg^dx) - J2 



x] 



<es,D 



(24) 
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where 



= Px/fT^x = ^xiX^v^) is the 4-velocity of packet x (x = a, 6) and, finally, q is the 
4-momentum of the virtual neutrino. Let us now define the symmetric overlap tensors 
for the source and detector: 

gfjf = T^^ and 3fJf = ^ T^^ 

xG5 x&D 

One can prove that these tensors are positive definite provided ^ Q (Vx); thus, there 
exist the positive-definite tensors J?^'^ and such that {^s,d)xu ~ the 
matrix form, 3?^^^ = ||3fJ^^'^|| = g'^'^g, where 3?^,^ = ||3fJ^';i|| and = \\g^,u\\- Obviously, 
\^s,d\ > and \^s,d\ = |3f^s,d|~^- These facts allow us to compute the integrals (1211) 
explicitly. For this aim, we define the quantities 

6s,,{K) = (47r)-2|3i,,,|-i/2exp (^-i^^^^iT^K,^ , (25) 
Kd = Kd E ^-^-^ = E - ' (26) 

©s,. = E (^-'^"^ - ^S' ^^^v) (27) 

K,X' 

(here x, x'eS, D and i^' is an arbitrary 4-momentum). With this notation, we find the 
following compact expression for the integral f l24l) : 

YsM) = (27r)''5s,d {qTqs,d) exp [-©s.d ± i {qTqs,d) ^s,d] ■ 

Let us clarify the physical meaning of functions ( !25l) - (!27j) . In the plane- wave limit 
(cTjf —7- 0) the factors 6s {q — qs) and 6d {q + qd) become usual 5 functions responsible for 
the exact energy-momentum conservation in the source and detector vertices, whereas 
at nonzero a^, they only lead to an approximate energy-momentum conservation: the 
probability of the process Is® Id Fs®Fd is strongly suppressed at a small disbalance 
in the 4-momenta of the interacting packets and the allowed disbalance is defined by 
the tensors S?^'' and 3?^*^, i.e., ultimately, by the momentum spreads of the packets. 

The functions exp (— ©s) and exp {—&d) are the geometric suppression factors 
conditioned by a partial overlap of the wave packets in the spacetime region of their 
interaction. This can be seen after converting (!27|) to the forml||J 

es,d = E (^^ - Xs,d)^ (x. - Xs,d), (28) 

X 

and taking into account that &s,d and Xs^d are invariants under the group of 
transformations 1 — > x° + r^^, cc^^ 1 — > + v^t,^ with arbitrary real parameters 



In this derivation we have used the translation invariance of the functions 6 



s,d- 
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r^. The latter symmetry allows, for each packet xgS", D having a nonzero velocity w^, 
the vector in fl28|) to be replaced with the vector 

b^= {x^- Xs,d) - [n^{x^- Xs,d)]n^ {n^ = v^/\v^\) (29) 

(whose absolute value, \b^\ = |n-^ x {x^ — X s,d)\, is the minimum distance between 
the classical world line of the center of the packet x and the point Xg^^j, and the zero 
component a;^ to be replaced with the time 6° of maximum approach of the packet 
center to the point Xg ^, this time being equal to 

€ = K - X%) - \v^\-'n^ (x^ - X,,,) . (30) 

The 4- vector = (6° , b^), built from ( l29i) and ( l30l) . is a relativistic analog of the usual 
impact parameter, so it is natural to call it the impact vector. The 4-vectors and X^ 
can be called, accordingly, the impact points. If ^ (Vx) then, after the substitution 
I — )■ b^, the expression (l28l) becomes 

= [(r^ - 1) + bi] = 

In the last equality, the contribution from each packet is written in its own rest frame (in 
which |b* I = I a;* — X*^|), and hence the temporal limitation v^^^ can be bypassed. 
The physical meaning of the factors exp (— ©s) and exp (— ©d) is now transparent: the 
interaction of in and out packets is unsuppressed (©s,d ^ 1) if all impact parameters 
|6* I are small relative to the effective dimensions of the packets k (~ Xjoy^. In the lab. 
frame, the suppression of "unlucky" configurations of the packets' world lines are defined 
by both the space and time components of the impact vectors 6^. The contribution of 
the time components 6° is less important for nonrelativistic packets (with ~ 1), but 
is very essential for ultrarelativistic packets (with ^ 1). 

Our consideration suggests that the impact points Xg and X^ identify the spacetime 
position of the regions of effective interaction of the packets in the source and detector, 
respectively. The more intensive interaction of the packets, the closer their world lines 
are placed with respect to the impact points. The condition that the interaction regions 
in the source and detector are macroscopically separated is equivalent to the macroscopic 
separation of the impact points. The world line configurations and impact point 
coordinates have no concern with the dynamics governed by the interaction Lagrangian, 
being uniquely defined by the initial (final) coordinates, group velocities and effective 
dimensions of the asymptotically free in (out) packets. But, of course, the full probability 
of the process Is®Id — ^ Fg^F^ is defined by both the suppression factors exp (— ©s) and 
exp {—&d) and interaction dynamics. 

It is evident now that the spacetime remoteness of the initial 4-coordinates of the 
in and out wave packets from the interaction region means simply that 

< \xa - Xg^dl > max„ (a^^) , aElg^d, 

xl > X^a, \xh - Xg^dl > maxfc (a^^) , beFg^d- 

Therewith, for the applicability of the CRGP model it is necessary that the components 
X^^ — x'^ remain finite in magnitude and satisfy the conditions ( |20l) . It is meaningful to 
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note that the requirement of spatial remoteness is in fact unnecessary for specification 
of the asymptotically free initial and final states. Indeed, some of the packets (e.g. a 
target nucleus or decaying meson in the source vertex) can be at rest in the lab. frame 
and thus they must be spatially close to the corresponding impact point; otherwise, 
the resultant interaction amplitude will be strongly suppressed by the geometric factor 
exp (-©s) or exp (-©d)- 

4. Calculation of a macroscopic amplitude 

As a practically important example, we consider the charged-current induced production 
of the charged leptons £^ and £^ (£a,/3 = e, /i, r) in the process 

J,©/rf^F; + £+©F^ + £-. (31) 

We assume for definiteness that all the external substates Is, Id, F's and F'^ consist 
exclusively of (asymptotically free) hadronic wave packets. Consequently, ii a ^ (3, the 
process (13T]) violates the lepton numbers and Lp that is only possible via exchange 
of massive neutrinos (no matter whether they are Dirac or Majorana particles). In 
the lowest nonvanishing order in electroweak interactions, the process fl3T|) is described 
by the sum of diagrams of figure |H Let and be the impact points defined by 




Figure 4. A macroscopic Feynman diagram describing the process (PT|) . 

eq. f l2B]) . We require these points and by that the effective regions of interactions in the 
source and detector (the areas symbolically outlined in figure H] by dashed curves) to be 
macroscopically separated and the conditions 

^a«^M (Va €/,,,) and x° » X", (V6 G F,,,) 

to be fulfilled. Hence the incoming and outgoing states are thought to be direct products 
of free one-packet states ([2]), each normalized according to the relation (ITOl) . 

In the framework of the standard model (SM) phenomenologically extended by 
inclusion of a neutrino mass term, the quark-lepton blocks of the diagram of figure H] 
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are described by the Lagrangian 

Cwix) = [je{x)W{x)+jg{x)W{x) + H.c] , 

where g is the SU{2) gauge couphng constant, and jq are the lepton and quark weak 
charged currents, 

= Y,v:M^)on^{x), j!;{x) = 5^\/';,g(x)0V(x), 

ai qq' 

Vai (a = e,/i, r, i = 1,2,3) and V^^, (g = u,c,t, q' = d,s,b) are the elements of the 
neutrino and quark mixing matrices {V and V' , respectively), £q,(x) and i/six) are the 
charged lepton fields, and O'^ = ■y'^ {1 — 75). The customary notation is used for other 
fields and Dirac 7-matrices. The normalized dimensionless amplitude of the process (15T|) 

(out|S|in) ((in|in) (out|out))^"'"^^ =^ A/sa 
is given by the fourth order of the perturbation theory in the coupling constant g: 

Ape. = ^ {^^^ {Fs®Fa\T j dxdx'dydy' : 3,{x)W{x) : : 3q{x')W\x') : 

X : 3\{y)W\y) : : 3l{y')W{y') : §,|/.©/,). (32) 

Here S/j = exp \i f dzC^iz)] , Ch{z) is the Lagrangian of strong and electromagnetic 
interactions responsible for nonperturbative processes of fragmentation and hadroniza- 
tion; T and : . . . : are the standard symbols for the chronological and normal ordering 
of local operators. The normalization factor A/" in the CRGP approximation is given by 

AT^ = (in|in)(out|out) = JJ 2E^V^{pJ. (33) 

The assumed narrowness of the wave packets in the momentum space, macroscopic 
remoteness of the interaction regions in the source and detector, and the consideration 
of translation invariance allow us to represent the hadronic part of the amplitude ([32]) 
in the factorized form 

(F;©F^|r [: ji;{x) ■. ■. 3l''{y) ■. §.] IL®!,) = jnPs)Jr{pD)ii\ 

X [ H e-'P^^^MPa, Xa -y)][ll e'^'^'^rbiPb, Xb - y) , 

where Js{ps) and Jdipo) are the c-number hadronic currents in which the strong 
interactions are taken into account nonperturbatively, and ps and pr, denote the sets 
of the momentum and spin variables of the hadronic states. Now, by applying Wick's 
theorem and the known properties of the leptonic wave packets, the amplitude fl52]) can 
be rewritten in the following way: 

^^- = ^11 ^P^^rn{Pp)0^'Ku',', ({p., 0''v{p,)J^V:^, (34) 
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where we have introduced the tensor function 

Gi.V.({p.,xJ) = I ^V,(g)A..,(g-j,^)A^(g)A,,,(g + p„)V.(g). (35) 

Here Vs(g) and Vd(g) are the overlap integrals discussed in detail in section [3l 

A"'(g) = i {q ~ nij + iO)~^ 

and A^y are the propagators of, respectively, the massive neutrino Uj and W boson (the 
explicit form of A^^ is not used below), mj and m]y are their masses, v{p^) and u{Pis) 
are the Dirac bispinors describing the leptons i'^ and i^, respectively. Here and below, 
the spin indices and arguments of the functions J's^d are omitted for short. 

Let us examine the amplitude ([31]) at large values of \Xs — X^l- The integral 
(!35|) can be evaluated by means of the theorem of Grimus-Stockinger (GS) [9]. Let 
$ = $(g) be a thrice continuously different iable function such that $ itself and its first 
and second derivatives decrease not slowly than as \q\ — )■ oo. Then, according to 
the GS theorem, in the asymptotic limit of L = |L| — )■ cxd, 

dq ^{q)e"i^ f (^L/L) exp (^v^L) + O (L-3/2) at s > 0, 

(27r)3 s-q^ + iO^ \ ^ O {L~^) at s < 0. 



The integrand in ( !35ll satisfies the formulated requirements. The integral over go? which 
remains after applying the GS theorem, can be evaluated by the regular saddle-point 
method. Here we describe only the result of this lengthy calculation performed in 
the ultrarelativistic approximation, which means that g° ~ — g° ^ m^, j = 1,2,3. 
The stationary saddle point go = Ej (which can be naturally treated as the effective 
energy of the virtual neutrino uj) is found as a series in powers of the small parameter 
r, = m]/{2ElM 

Yl Yl 

E,=E4l-nr, + 0{r])], = ^=yI^ 

Y'^ = W/qs, - ^^qdu, R = + K) ^m^- 

The effective momentum and velocity of the virtual neutrino are then found as Pj = Pjl 
and Vj = Pj/Ej = Vjl, respectively, where 



P, = Je] - m] = i^;, [1 - (n + 1) r, + O (r|) ] 



and 



r, + O (r|) . 

As is easy to see, Ej = Pj = E^ = q^ = —q^ in the limit of mj = and assuming 
the exact energy-momentum conservation. But, in the general case, the effective 4- 
momentum pj = {Ej,Pj) is determined by the mean momenta and effective dimensions 

% We limited ourselves to the first order of the perturbation expansion in Vj. However, the next-order 
corrections are, in fact, needed to define properly the range of applicability of the obtained result. 
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of the external wave packets involved in the process 0311] . Finally, by introducing the 
notation 

n,(T, L) = I {E,T - P,L) + ?S (p^.T - E,Lf , (36) 

3 

Q = Xsqs + Xdqd, L=\Xa-X,\, T = X° - (37) 
2^2 r\ -V2 



(38) 



2),- = D,. 1 + — ^ ^ , 2), = 

' ' V Pf J ' V2R 

we arrive at the saddle-point estimate of the function f l35l) : 

^.g-n,(T,L)-ie 

This formula can be, and must be, somewhat simplified by putting rj = everywhere 
wherever it is not a factor multiplying L or T (values which can be arbitrary large). 
Then the 4-vector pj is replaced by pu = {Ey^p^^) = EJ. Taking into account the 
aforesaid, the complex phase (136!) can be written in an explicitly invariant form 

n,{T, L) = t{p,X) + [{p,Xf - m]X^] , X = X, - X,. (40) 

Comparing the factor ip* = e^^^ in (!39|l with the generic CRGP wavefunction (!2T|) . we 
conclude that ip* can be treated as the (outgoing) neutrino wave packet in which the 
role of the parameter a is played by the function = \/2'Dj/Tj (Tj = Ey/mj). Since 
the latter is a complex-valued function, the neutrino wave packet spreads with increase 
of L = \X\. The spread effect appearing at very large L and leading to both the overall 
suppression of the amplitude (IMIl and modification of the relative contributions with 
different j into can be of definite interest for neutrino astrophysics. However, 
in this paper, we limit ourselves to the analysis of the probability of the process 
fl3T|) under "terrestrial" conditions, for which it is a fortiori possible to judge that 
E,jL < {VjEj/2Tijf and, consequently, to put ~ Dj ~ l/y/2R = D and 

~ V2^/Tj = l/iTjy/R). Hence, the relative energy-momentum uncertainty of 
the ultrarelativistic neutrino packet, SEj/Ej ~ 6Pj/Pj ~ D/E^, is small and is 
determined by the momentum spreads of the external in and out packets. Of course, 
the mean position of the neutrino wave packet evolves along the "classical trajectory" 
L = VjT, the quantum deviations from which, 6L, are suppressed by the factor 
exp{-22)2 [5LyT] + {l5Lf]}. 

Now, by applying the identity P-p^P^ = P-U^{py)u^{p^)P+, in which u_(pj,) is 
the Dirac bispinor for the left-handed massless neutrino and P± = |(1 ± 75), we define 
the matrix elements 

g2 

Ms = ^u4p,)Js^A^^,{p, + p^)0^'u{p^), 

m: = ^^v{p^)o^'A,,,{p^-pp)r/u^{p,), 

describing, respectively, the production of a real massless neutrino v in the reaction 
Is — )• F'sC-t^ and its absorption in the reaction ul^ — )■ F'^f-^. Then, taking into account 
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the above-mentioned results, we obtain the final expression for the amplitude 

^(27r)3/WL 



As is evident from the derivation of this formula and its structure, it is valid for 
essentially any class of macrodiagrams with exchange of virtual massive neutrinos 
between the source and detector vertices, unless we do not specify the explicit form 
of the matrix elements Mg and M^i. To obtain similar result for the macrodiagrams 
with exchange of virtual antineutrinos, one must replace (besides the matrix elements) 
V by Ft (i.e., ^ K„ V^, ^ V^*^). 

5. Probability and count rate. 

It can be shown that 

\^sAPu)\' = i27fY5s,d{PuTqs,d)ys4, (42) 

where 6s^d (the "smeared" 6 functions analogous to the functions 6s^d) and Vsa (the 
effective four-dimensional overlap volumes of the external packets) are given by 

5.,4K) = (27r)-2|SR,,,|-V2 exp (^-1 ^'^^.K.K,^ , (43) 
V.,rf = fdxH \^Ap^,x^-x)\\ (44) 
Thus, from (14T|) we obtain the microscopic probability of the process (13T]) . 



X 



(27r)3L2 



dependent on the parameters a,^, coordinates and mean momenta of all packets 
participated in the reaction. The probability (H5|l is small if the product of the overlap 
volumes (Hill . 

V,V, = (7r/2)^ il^Mdir'^' exp [-2 (6, + &,)] , 

vanishes, i.e. if the external wave packets in the source and detector do not overlap in 
the spacetime regions surrounding the impact points Xg and Xa- 

Let us note that the 4-vector p^, is a function of and cr^, and Pu = = —Id in 
the plane- wave limit {a^ = 0, Vx). Thus, for sufficiently small a^, 

5s{P. - qs)5d{Pu + Qd) ^ 5s(0)5d(0) = (27r)-^ il^Mdir'^' . 

But what controls the approximate equality of the 4-momentum transfers and — g^? 
To answer this question, it is useful to rewrite (HSj) in the fashion used by Cardall for 
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a similar purpose [TT]. By using the explicit form of the functions Sg^d and D, one can 
prove the following approximate relation: 

2V^^6s {p. -qs)6d {p. + qd) F{p,) = J dE% (p'. - Qs) ^d (p. + qd) F{pl) , (46) 

which is valid with the same accuracy with which the amplitude f l4T|) itself has been 
deduced (that is, with the accuracy of the saddle-point method). Here F{py) is an 
arbitrary slowly varying function and = {El,pl) = E'J. With the help of the 
relation (l46l) the microscopic probability (l45l) transforms to 

(27r)^(5,(p.-g,)V,|M,p {27c)%{p, + q,)V,\M,\^ 

^ 'j:v:,V,,e-^^^^^'^\ (47) 



2v/^(27r)3L2 

where we have omitted the prime on the integration variable E^,, but now E^ (as well 
as = E^l) is in no way related to the parameters of the external wave packets. 
Expressions (H5l) and (H71) are equivalent within the adopted accuracy, but from (H71) it is 
apparent that the energy- momentum conservation is governed by the factors Ss{pu — qs) 
and 6d{Pu + qd) which, at sufficiently small a^, could be substituted by the usual S 
functions. 

The probability (H7|) is the most general result of this paper. However, it is too 
general to be directly applied to the contemporary neutrino oscillation experiments. 
In order to obtain the actually observable quantities, the probability P7|) should 
be properly averaged/integrated over all the unmeasurable or unused variables of 
incoming/outgoing wave-packet states. Let us call this procedure the macroscopic 
averaging. Clearly, such a procedure can only be realized by taking into account the 
physical conditions of the real experimental environment. For these reasons and in this 
sense, further analysis becomes model-dependent. 

As a simple but realistic example, we consider a thought experiment in which it 
is assumed that the statistical distributions of the incoming packets a G Is^ over the 
mean momenta, spin projections and spacetime coordinates in the source and detector 
"devices" can be described by the one-particle distribution functions /a(Pa, Sa, Xa). It is 
convenient to normalize each function fa to the total number, iVa(a;°), of the packets a 



at a time 



Sa 



For clarity purposes, we must define (or rather redefine) the terms "source" and 
"detector" which were so far used for designating the blocks of the macrodiagrams. In 
what follows we will use these terms and notation S and V both for the corresponding 
devices and, more abstractly, for the supports of the products of the distribution 
functions /„ in the spacetime variables (namely, S = suppj^^j Yla fa, ciElg and similarly 
for V), which are assumed to be finite and mutually disjoint within the space domain. 
We further suppose that the effective spatial dimensions of S and V are small compared 
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to the mean distance between them but very large compared to the effective dimensions 
(~ cr~^) of all wave packets moving inside S and V. For definiteness, we also approve 
that the experiment measures only the momenta of the secondaries in V, and (owing to 
the large distance between S and V) the background events caused by the secondaries 
falling into V from S can be neglected. Lastly, we accept the detection efficiency to 
be 100%, though the formalism allows a straightforward account for the real efficiency, 
acceptance, etc. With all these assumptions, the macroscopically averaged probability 
( H7D reads 

f TT dXgdpJajp^, Sg, Xg) r -j-r dXb[dp^] 

J 11 (2vr)32E.V. J 11 (2.)32i^,V. ' 

X J dE,{2TT)*6s{pu-qs)\M,\\27r)%{p, + qd)\Md\^ 
2) 



X 



20F(27r)3L2 



J2^:,V,,e-^^^^^'^ . (49) 



Here and below the symbol Xlspins denotes the averaging over the spin projections of 
the initial states and summation over the spin projections of the final states in S and 
D. The symbol [dpf^] indicates that integration in the variable pj, is not performed, i.e. 
J[dph] = dpf,. With regard to the normalization conditions (l48l) . it is easy to recognise 
that fj49|) represents the total number, dN^p, of events recorded in V and consisted of 
the secondary wave packets h^F^ having the mean momenta between and + dp^^. 
Under additional assumptions, the somewhat unwieldy expression ( H9i) can be simplified 
in a few steps. 

An approximate multidimensional integration over the spatial variables in (H9ll 
can be performed using the integral representation for the overlap volumes (jHj) 
and taking into account that the distribution functions /a, as well as the factors 
^-vij{T,L)-n*{T,L) I -£2 assumed to vary at large (macroscopic) scales, whereas the 
integrand fix l'^^^ (Pjf -^^^ ~ (l44l) is essentially different from zero only if the 

classical word lines of all packets x pass through a small (not necessarily microscopic) 
vicinity of the integration variable (let the latter be x for and y for V^). Hence, 
neglecting the edge effects, it is safe to replace x^ by x {y) for x G 5 {D) in the 
mentioned slowly varying factors. Then, as is seen from (l26l) . = x and = y. The 
remaining integrals over the variables yield the factor rij^ese-D which cancels the 
same factor in the denominator of the integrand in ( l49l) . As a result, we can rewrite 
(|19|) as follows: 

2 



spins 



2) 



dN^p = J2 I dx I dy I d^is I d^d / dE, l67rV2|^ _ ' (^0) 
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where we have defined the differential forms 



dPafa{Pa,Sa,x) 



n 



dPb 



n dPafa{Pa,Sa,y) -j-r 



The phase Qj (T, L) in 



i2n)%{p, + qd)\Md\ 



(51) 
(52) 



501) is still defined by equation ( l36ll or (HOl) . in which, now 



T = Xq = Uq — Xq and L = \X\ = \y — x\. It must be underlined that the considerations 
used for deriving ( ISOl) were based on the additional restrictions, which may not be fully 
adequate to the particular experimental conditions^ Thus the comparatively simple 
but approximate formula ( !50|) is not entirely equivalent to the more sophisticated and 
general result ( H9ll . 

For further simplification, we perform in f l50|) integration in time variables and 
. This integration can be easily done in assumption that, during the experiment, the 
distribution functions fa in S and T> vary slowly enough with time so that they can be 
modelled by the "rectangular ledges" 

/a(Pa, Sa] x) = 9 (x° - X?) 9 {x9. - X°) Ja{Pa, ^a] x) for aG/„ 



fa{Pa, Sa, V) = 9 (y° - y?) 9 (y° - 7/°) /„(p„, Sa, V) for OG/rf. 



(53) 



In case of detector, the step functions in (!53|l can be thought as the "hardware" or 
"software" trigger conditions. The periods of stationarity = X2 — x^ and Td = 1/2 — Hi 
can be astronomically long, as it is for the solar and atmospheric neutrino experiments 
{ts ^ Td in these cases), or very short, like in the experiments with short-pulsed 
accelerator beams (when usually < r^), but we anyhow assume that the time 
intervals needed to switch on and switch off the source (detector) are negligibly small 
in comparison with Tg (r^). Within the model (l53il . the only time-dependent factor in 
the integrand of fl50|) is e~^^~^i . So the problem is reduced to the simple integral 



dy' 



^-nj{yO-xO,L)-nUyO-xO,L) 



2^ 



Td exp [iipij - Si_ 



In this relation, we have adopted the following notation: 
exp (— 



S,. 



i,i'=i 



lerf 



2D xi- - yj; + 



L 



(v, - v,)T)L 



where 



27rL 



27tT)L 

EyLij 



AE^ 
4D 



vrn 



2DL,-,- 



(54) 

(55) 
(56) 



1 

Vi 



+ Moreover, we implicitly used the current experimental constraints on the neutrino masses, which 
suggest that lm{V,j + fl*) noticeably vary on the macroscopically large scales cx £'y/|m| — and 
Re(rij +^*) - on the scales much larger than L^. We have to remind, however, that such conclusions 
are based on the QM (rather than QFT) analyses of the existing data. 
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Ami = - ^^^J = ^^- 

/"^ 1 1 2 

Ierf(2:) = / dz'eri{z') + —= = zeri{z) + —=e , 

and erf(2) is the error function. For a more realistic description of the accelerator beam 
pulse experiments, the model ( !53l) could be readily extended by inclusion of a series of 
rectangular ledges followed by pauses during which fa = 0. We will, however, proceed 
with the simplest case, which reproduces the most significant effects. Then substituting 
(l5ll into ([50]) we obtain: 

dN., = r,Y:JdxJdyJd^J d^, J dEj^^^^^ (57) 

= j dx j dy j d^, j da.vVafiiE,, \y - x\). (58) 

The differential forms d^s^d in fl57|) are defined according to equations f lSTj) and fl52l) 
in which the distribution functions should be substituted by Equation fl58|) is 
written by applying the identity 

E_^^£^^d^^ _ d^^sdp^ d^ 

. 4(27r)3|i/-a;|2 ~ ^ (27r)^2E,\y - x\^dn, ^ ^2E, 

spins spins S 6 spins G 

where Vg and Vx> are the spatial volumes of the source and detector, respectively. The 
differential form d^^, is defined in such a way that the integral 

Vs dE.-'''' ^^J 2(2vr)3|^-xP ^''^ 

spins g A \ / ,^ 

is nothing else than the flux density of neutrinos in V, produced through the processes 
Is — )■ F^i'^u in S. More precisely, it is the number of neutrinos appearing per unit time 
and unit neutrino energy in an elementary volume dx around the point x E S, travelling 
within the solid angle dfly about the flow direction I = [y — x)/\y — x\ and crossing a 
unit area, placed around the point y eT> and normal to I. The quantity da^x> is defined 
in such a way that 

spins G U 

represents the differential cross section of the neutrino scattering off the detector as a 
whole. In the particular (and the most basically important) case of neutrino scattering 
in the reaction va — ?■ -F^^^, provided that the momentum distribution of the target 
scatterers a is sufficiently narrow, the differential form da^xi becomes exactly the 
elementary differential cross section of this reaction multiplied by the total number 
of the particles a in V. 

Now let us address the last subintegral multiplier of ( 158|) . given by 

Val3{Eu, L) = ^ VaiVajV^jVi3iSij CXp (iLfij - £/^j) . (61) 
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This factor coincides with the well-known QM expression for the neutrino flavor 
transition probability, provided that Sij = 1 and = 0, and thus it can be considered 
as a QFT refinement of the QM result. However a probabilistic interpretation of the 
function Vap can be only provisionally true, because the factors S^j and s^ij in (16T|) 
involve the functions S) and n strongly dependent on the neutrino energy E,y and external 
momenta p^; all these (except for the momenta of secondaries in V) are variables of 
integration in (|58|) . As a result, the factor Va/s, as function of a and (3, does not satisfy 
the unitarity relations 



which are a commonplace in the QM theory. The point is that the domains and shapes 
of the functions D and n are essentially different for each of the nine leptonic pairs 
{C-aAp)- These differences are governed by the kinematics of the subprocesses in S and 
V (in particular, their thresholds), that is, eventually, by the leptonic masses (mg, m^, 
nir) and by the momentum spreads (cTg, a^, CFr) of the leptonic wave packets, which are 
not necessarily equal to each other, perhaps even within an order of magnitude. The 
probabilistic treatment of Vap is even more problematic in real-life experiments, because 
the detector event rate (with £^ appearance in our case) is defined by many subprocesses 
of different types in the source and detector. For example, in the atmospheric and 
accelerator neutrino experiments, the major processes of neutrino production are in- 
flight decays of light mesons (vr^25 K^2^ K^^, K^s, etc) and muons, and neutrino 
interactions with detector consist of an incoherent superposition of exclusive reactions 
of many types, - from (quasi)elastic to deep inelastic. 

One more technical drawback to interpretation of (!6T]) is the dependence of the 
function Sij (which will be referred to as a decoherence factor) on the four "instrumental" 
time parameters x^, x^, Vi and 7/2 . So far we have made no assumption concerning a 
"synchronization" of the time windows {xi,X2) and (yi,l/2)- Thus, it is no wonder that 
the decoherence factor turns to be vanishingly small in magnitude if these windows 
are not adjusted to account that the representative time of ultrarelativistic neutrino 
propagation from 5 to P is equal to the mean distance, L, between S and V. 
Before discussing the role of the decoherence factor, we perform one more, and the 
last, simplification of the formula for dNap-, again using the requirement that the 
characteristic dimensions of S and T> are small compared to L. Under certain conditions, 
this allows us to replace \y — x\ in the integrand of by L and the differential forms 
d^y and da^-o by their averages, d^^ and (ia^x), over the spatial volumes of S and T>, 
respectively. Finally, we arrive at the approximate formula 



Its range of applicability is in general much more limited than that of fl58|) . as a 
consequence of additional restrictions implicitly imposed on the distribution functions 
/q, absolute dimensions and geometry of S and V. These issues are not discussed 




(62) 
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in this paper but must be the subject of special attention in the neutrino oscillation 
experiments. 

Let us now return to the decoherence factor, limiting ourselves to a consideration 
of "synchronized" measurements, in which = -FTs/2, 1/12 = With certain 

technical simplifications, the factor ( l55ll can be expressed through a real- valued function 
S{t, t', b) of three dimensionless variables, namely 

2t'S{t, t', b) = exp (-62) Re [lerf {t + t' + ib) - lerf [t-f + ib)] . 

Figure |5] shows the numerically evaluated 3D plot and 2D density plot of the function 
So{t,t') = S(t,t',0). It can be proved that < So(t,t') < 1 for any t,t' > and 




Figure 5. A 3D plot (left) and a 2D density plot (right) of the decoherence function 
So{t, t'). The darker regions in the 2D plot correspond to the smaller values of SQ{t, t'). 



So{t,t') < t/t' for t' > t. This implies that the mean count rate in the detector, 
dRajs = dNais/Td, dccreases with a rise of the ratio r^/ > 1. The reason is apparent: the 
number of events recorded in V cannot be larger then the number of neutrinos emitted 
from S. A less obvious conclusion follows from the inequality So{t + 6t,t) > erf (St) 
valid for 6t > 0: the count rate is not suppressed at sufficiently large Ts/rd- Clearly, this 
condition is over-fulfilled in the solar and atmospheric neutrino experiments but may be 
violated in the accelerator and perhaps reactor experiments. Two dashed lines in the 
right panel of figure [5] separate the regions in which Sq < 0.5 (r^ < 2Td) and > 0.995 

{Ts > Td + 2/D). 

In the particular case t = t', which is of interest for accelerator experiments, the 
function 5'o(t,t) converges to unity only at very large t (in practice for t > 100). Thus, 
in order to set So(t,t') = 1, the detector exposition time must be either sufficiently 
small in comparison with or (at ~ Tg) much larger than the characteristic time 
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scale Tj, = 1/ min(S}), where the minimum should be taken over the phase subspaces of 
the process fl3Tl) . which are responsible for a significant contribution into the count rate. 

On the other hand, the observed strong dependence of the common suppression 
factor So(t,t') on its arguments at t < t' provides a potential possibility of an 
experimental estimation of the function D (or, rather, of its mean values within the 
above-mentioned phase subspaces), based on the measuring the count rate dRai3 as a 
function of and (at fixed L), and comparing the data with the results of Monte- 
Carlo simulations. The optimal strategy of such an experiment should be a subject of 
a dedicated analysis. 

The profiles of the function S{t,t',b) calculated at various fixed values of b are 
displayed in figure [61 It appears that the shape of S{t,t',b) as a function of t and t' 
becomes more and more complicated with an increase of b. For 6 > 3 — 4, S{t, t', b) 
rapidly oscillates around zero, highly suppressing the "non-diagonal" (with i ^ j) terms 
in the factor Va/3iE^, L). Figure [7] shows the dependence of the function S'(t, t, b) on the 
parameter 6, evaluated at various fixed values of t. At large t, this dependence has a 
quasi-periodic character superimposed by a fast decrease of S{t, t, b) with increasing 
b. At very large t, the function S{t,t,b) becomes nearly independent of t, slowly 
approaching the asymptotic behavior 

5(t,t',6) ~ exp(-6^) {t,t' 

In this asymptotic regime, the "probability" (|6Tll takes on the form already known from 
the literature (see, e.g., |6l[13] and references therein), 

V.^iE,, L) = J2 K^Va,V;^Vp, exp [i^,, - - ) , (63) 

but with the essential difference that the factors s^ij and do depend (through the 
functions 25 and n) on the neutrino energy and momenta of the external wave packets. 
This dependence drastically affects the magnitude and shape of these factors if at 
least some of the wave packets have relativistic momenta (that is always the case in 
the contemporary neutrino oscillation experiments). As follows from our analysis, for 
sufficiently small and/or hierarchically different momentum spreads a^, the functions 
s^ij and ^iSij may vary in many orders of magnitude through their multidimensional 
domain. 

The factors exp {—■^^) (with iy^j) suppress the interference terms in f l63|) at the 
distances exceeding the "coherence length" 

when the neutrino wave packets ip* and ip* are strongly separated in space (due to the 
difference in their group velocities) and no longer interfere. Clearly L™*^ — cxd in the 
plane- wave limit. 

The suppression factors exp (— (i^j) work in the opposite situation when the 
external packets in 5 or P (or in both S and V) are strongly delocalized (in the plane- 
wave limit - uniformly distributed over the whole space). The gross dimension of the 




Figure 6. Profiles of tfie decolierence function S{t, t' , h) calculated at nine values of b 
shown on the panels. 

neutrino production and absorption regions in S and V is of the order of The 
interference terms vanish if this scale is large compared to the "interference length" 

Ljf = 1/(4AE,,) = 2L,,/(7rn). 

In other words, the QFT approach predicts vanishing of neutrino oscillations in the 
plane-wave limit. In this limit, the flavour transition probability does not depend on L, 
Ejj and neutrino masses and becomes |V^iP|V^ip. Thereby, a nontrivial interference 
of the diagrams with the intermediate neutrinos of different masses is only possible if 
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Figure 7. The function S{t,t,b) versus b at fixed values of t shown in legend. The 
asymptotics of S(t, 6) at t — > oo and finite b is shown for comparison. 



D ^ 0. Our analysis of the generic subprocesses 1 — )■ 2, 1 — )■ 3 and 2 — )■ 2 shows that 
2) 7^ if in both vertices of the macro diagram describing the process fl3Tl) there are 
at least two interacting wave packets x (no matter in or out) with cr^ ^ 0. The same 
requirement unavoidably leads to the vanishing of the non-diagonal terms, when the 
mean distance between S and V becomes large enough in comparison with the coherence 
lengths L^°^. As a result, the range of applicability of the standard QM formula for the 
neutrino oscillations probability is limited by rather restrictive conditions, 

in which the angle brackets symbolize an averaging over the phase subspace of the 
process fl3T|) which provides the main contribution into the measured count rate. We 
recall that the conditions ( l64l) . as well as the relation f l63l) . were obtained under a number 
of assumptions and simplifications, which are not necessarily adequate to fully represent 
the real-life experimental conditions. Our consideration suggests that in the analysis and 
interpretation of real data, one should take into account not only the conditions flM|) . 
but also the operating times of the source and detector, their geometry and dimensions, 
explicit form of the distribution functions of in packets and other technical details. 
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6. Conclusions 

In this paper we have considered a covariant theory of wave packets in which the free 
field states are constructed as superpositions of one-particle Fock states and are reduced 
to the Fock states in the plane- wave limit. The mean momentum and effective spatial 
volume of such a wave packet are exact integrals of motion, in spite of the natural 
dispersion of the packet with time. The mean position of the packet evolves along the 
classical trajectory. As a simplest model, satisfying all requirements of the formalism, 
we have studied relativistic Gaussian packets (RGP). The conditions are found which 
allow one to neglect the spreading of RGP (CRGP approximation) and use the latter 
for description of the asymptotically free particle states in the S'-matrix formalism of 
QFT. 

The formalism has been implemented for the calculation of Feynman diagrams 
with massive neutrino exchange between two macroscopically separated vertices and 
with the wave packets as external legs. As a typical generic example, we have studied 
the amplitude of a lepton flavor violating process or rather a class of the processes 
( I3TI) with production of the two charged leptons ia and ijs in the distant parts of the 
macrodiagram. The amplitude is proportional to the sum of the products of the matrix 
elements describing the subprocesses of neutrino creation and absorption and the factor 
ip* I L, which can be treated as a spherical wave of neutrino Vj with a mass my The factor 
turns to be of the same form as the generic CRGP wavefunction, but its "dispersion" 
is a relativistic-invariant function of the neutrino energy Ey and momenta of external 
wave packets p^. As a common multiplier, the amplitude contains a suppression factor, 
arising due to the incomplete spacetime overlap of the external wave packets and, in 
addition, two non-singular factors, responsible for an approximate energy-momentum 
conservation in the vertices of the macrodiagram. 

Experimental detection of the flavor- violating process ( 1311) would provide a firm 
indication of the incompleteness of the standard model, and a detailed measurement 
of its probability would be an important source of information on the neutrino mixing 
parameters and squared-mass splittings. A Lorentz-invariant formula for the probability 
of the process ( 13T|) has been obtained. A macroscopic averaging of this probability 
leads to an experimentally measurable quantity - the differential number of events 
represented, after several simplifications, by a multidimensional integral of the product 
of three factors, see ( 1581) . The factors d^y and davx> are respectively related, but not 
identical, to the differential flux density (energy spectrum) of massless neutrinos from 
the source and to the differential cross section of the neutrino-detector interactions. 
The function Vafi is answerable for the neutrino flavor transitions. But in general 
this quantity does not have the properties of a probability, because it involves the 
decoherence factors dependent on the phase space variables Ey and ip,^ and on the 
masses and momentum spreads of the external wave packets, including the packets of 
the leptons ia and i^. 

Within a simple model for the distribution functions of in packets in the source and 
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detector, the magnitude of the common suppression factor 5*0 in Vap is determined by the 
operating times of the source (r^) and detector (r^), and by the spacetime width of the 
neutrino wave packet ~ 1/2), where 25 is a function of Ei, and p^. The suppression is 
small if Tg ^ Td or Tg ^ t^. In the opposite case, the strong dependence of Sq on the 
parameters Ts and provides a potential possibility of measuring an average value of the 
function 2) in a special-purpose accelerator experiment, which allows us to variate these 
parameters (or at least one of them). Such a measurement, even if rather rough, would 
be useful in design and data-handling of future precision experiments with neutrino 
factories, /3-beams and super-beams. The non-diagonal decoherence factors Sij ii^j) 
show in general a very involved behavior. In an asymptotic regime, Tg ^ r^, these 
factors cease to depend on and but continue to strongly suppress the "oscillation" 
terms in Vap- The conditions have been found at which the decoherence effects are 
small and the standard QM formula for the neutrino flavour transition probability is 
applicable. 
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